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Based on the generalized non-Markovian equations obtained earlier for a nonequilibrium one-particle distri- 
bution function and potential part of the averaged enthalpy density [Markiv B.B., Omelyan I.P, Tokarchuk M.V., 
Condens. Matter Phys., 2010, 13, 23005] a spectrum of collective excitations is investigated, where the potential 
of interaction between particles is presented as a sum of the potential of hard spheres and a certain long-range 
potential. 
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1. Introduction 

A number of investigations were devoted to the problem of constructing a consistent descrip- 

tion of kinetic and hydrodynamic processes in dense gases, liquids, and plasma. For instance, the im- 
portance of taking into account the kinetic processes connected with irreversible collision processes at 
the scale of short-ranged interparticle interactions was pointed out in |11]. Short-wavelength collective 
modes in liquids were investigated therein based on the linearized kinetic equation of the revised Enskog 
theory for the hard spheres model. 

In this paper we investigate a spectrum of collective excitations within a consistent description of 
kinetic and hydrodynamic processes in a system in which the potential of interaction between particles 
consists of two parts: the hard spheres potential and a long-range part. 

2. Transport equations 

Using the ideas presented in papers fljd] the nonequilibrium statistical operator consistently describ- 
ing the kinetic and hydrodynamic processes for a system of classical interacting particles was obtained 
m by means of nonequilibrium statistical operator method. Using this operator, a set of kinetic 
equations for the nonequilibrium one-particle distribution function fiip', t) - (h-^ip))' and the potential 
part of the averaged enthalpy density /j?'(f) - (K^^)' was obtained in the case of weakly nonequilibrium 
processes: 

^ flip: t)+'— flip; t)^-'—nfo{p)c2{k) [ Ap'fi(p';t) + ianh{hp)hf(t) 
at " m m J «^ 

t t 

-Jdp'l e'^''-'^ipnn(hp,p';t,t')fi(f;t']dt'- j e'^''-'^q>„Hik;p;t,t')hf\W, (D 
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3 r 

— 4"*(f) = j dp'inhnihp')fi[p';t) 

t t 

-Jdp'l e'^''-'^q)hn{hp';t,t')fi{p';t')dt'- j e'^*'-'''q)hh{ht,t')hf{t')dt', (2) 

— OO — CXD 

where in„,,(fc;p) = <Aj(p)^™^>o<I'^,',(fc) and \D.hn[k;p) = /dp'<^l"'n_^(p')>o<I'^^(p',p) are the normal- 
ized static correlation functions. 

cpnn(k;p,p';t,t') = J dp" {In(k;p)TQ{t,t')In(-hp"))Q<i]^{p" ,p'), 

(Ph„[k;p;t,t') = j dp'{I'^'[k)To[t,t')In{-k;p'))o^lhp',p), 

(P„hik;p;t,t') = </„(fc;p)ro(f,f')/f (-^)>o$^i(fc), 
(Phh{t,t,t') = </f (fc)ro(r,r')/f (-fc)>o$^i(fc) (3) 

are the generalized transport kernels (memory functions) describing kinetic and hydrodynamic pro- 
cesses. In l^^, the inner structure of generalized transport kernels for a consistent description of ki- 
netic and hydrodynamic processes was analyzed in detail. It was shown that they are expressed in 
terms of time correlation functions related to the basic set of dynamical variables, phase density n-j^ip) 
and potential part of the enthalpy density Iz?' along with the transport kernels describing diffusive 

and visco-thermal processes. Here, /i^(p) - f dre~^'"^hi (f, p) are the Fourier-components of microscopic 
phase density of particles number, hi(r,p) - T.'^L^Sip- pi)5{f- f/), fe"' - '^^^i^^-T^oS~^{k)h^ are 

the Fourier-components of the potential part of the enthalpy density, - |Zj^j=j'l'(l?/jl)e~^*^'^' and 

h-j^ - Xfli e~^^''' are the Fourier-components of the potential energy and particle number densities, re- 
spectively, k is the wave-vector. '^J^\(k) is the function inverse to the equilibrium correlation function 
9f,h(k] = <^!?'/z'"!)o, <...>o - f dTN---Qo(x^), where po is an equilibrium statistical operator. 7„(fc;p) - 
(1 - Po)iijv"fc(p) = (1 - Po)fiiip) and P^^k) = (1 - Pojiiw^l"' = d - Po)h'^^ are the generalized flows 
in linear approximation, To{t, t') - 6''"*^'"^"^"^'^^ tj^g evolution operator with regard to projection. 
Po is the linear approximation of the Mori projection operator constructed on the orthogonal dynamic 
variables hf' |8]: PqAj. = LiiAih^%)o%lihhf' + Lifdpfdp'{Ain_i^[p))o^^^^ It 

possesses the following properties: PqPo - Pq, Poil-Po) - 0, Pon-^ip] - nj.{p), Poh^-^^ - ^^^ip,p') 
is the function inverse to <l>j(p, p') - <nj(p)n_j(p')>o - n8{p- p')fo{p') + «^/o(p)/o(p')^2(fc)- It is equal 

to OZ^{p,p') - ^^f~^pt^ - C2(fc), where n- N/V, foip) - [fi!2nmf'^e~^^^ is the Maxwellian distribution, 
(i— l/fceT is an inverse temperature and fce is Boltzmann constant, czik) is the direct correlation func- 
tion related to the correlation function hzik): hzik) - C2(fc)[l - nczik)]'^. S(fc) = (nj_h_-f_)o denotes the 
static structure factor. It is important to note that dynamical variables h'^^ and n-j^ip) are orthogonal in 

the sense that (/i?*nj(p)>o - 0. 

Projecting the set of equations O onto the first moments of the nonequilibrium one-particle dis- 
tribution function = 1, ^a(p) = V2pa/2fcBr (where a = x,y,z), We[p] = VlPiip^ l2mk^T -2,12], 
one can obtain a set of equations for the averaged values of densities of particles number n-j^{t), momen- 
tum Ji{t), kinetic l^'^it) and potential /i™'(f) parts of enthalpy 1^, where the Fourier-components of 

fc k ^ 

the kinetic part of enthalpy density defined as li^" = e'?^ - {£^'^fi_-^)o{h-^h_-^)Q^h-^. For this purpose, we 
introduce the projection operator constructed on the eigenfunctions l^vip)) of the nonequilibrium one- 
particle function such that 3»|>P> =2:^^ J^v><^v|y>. Here, m'l^'v) = / dp^'[p)fo[p)^'y[p),wh\^e ^'vip) 
satisfies the conditions CV^l^v) - S^v and Xv \^v)Ci'v\ - 1- Then, let us act by the projection operator 
onto the set of equation l[2). Repeat this operation acting by the operator ^ - complementary 
to 3^. Then, substituting the unknown quantity from the second equation into the first one we obtain 
the necessary set of equations with separated contributions of kinetic and potential energies. Using the 
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Laplace transform, let us represent it in a matrix form: 

zaj^iz) - lG(k;z]a^(z] = -<%(f = 0))^ 
£g(^; z) is the matrix of memory kernels 

ZG(fc;2) = iOG(fc)-n(fc;z), 
where = [?ijf(z),J^(z),^?"(z),/zl"'(z)] is the column-vector. 



(4) 
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(6) 



are the frequency matrix and the matrix of transport kernels. The elements of the latter have the follow- 
ing structure: 



n^v (fc; z) = <^p|(^(fc; z) -H Z(fc; z)^ [z7 - ^£(fc; z)^] -i^£(fc; z)|^v>- 



(7) 



I denotes a unit matrix, (p(k;z) is the matrix whose elements are the generalized transport kernels 
(Pnn(k;p,p';z), (pfif,(k;p';z), (pfih(k;p;z), (pfifi{k;z) in the set of equations (T), l[2), and Z(fc;z) = iO(fc) - 
(p(k;z). Here, iQ(fc) is the matrix of static correlation functions iD.fihik;p), iO;;„(fc;p). For the sake of 
simplicity the dependence of (p{k; z), £(fc; z) on p, p' was omitted. As we can see from the structure of ele- 
ments of the matrices IOgC^) and n(fc; z), the contributions of kinetic and potential parts of enthalpy are 
separated. Herewith, a question arises regarding the study of time correlation functions and collective 
modes for Uquids based on the set of transport equations l|4). 



3. Spectrum of collective excitations 

Let us consider the system of kinetic equations (T), l[2} in the case where the potential of interaction is 
presented as follows: 

m\rij\)^':>'^H\Tij\) + <I>\\fij\), (8) 

where 't'^^i\fij\] is the hard sphere interaction potential, and O'dHjl) is the long-range potential. Taking 
into account the features of the hard sphere model dynamics |4] and the results of investigations |11, isj, 

one can separate Enskog-Boltzmann collision integral from the function (pnnik; p, p'; t, t') . Indeed, an 
infinitesimal time of a collision tq -i-O within an infinitesimal region cr+Aro, Aro ~ |To||p2-pi|/m ^ -i-O 
being a feature of the hard sphere model dynamics (a is the hard sphere diameter). Taking this into ac- 
count in the kinetic equation (1) we can obtain the kinetic equation of the revised Enskog theory for 
the hard sphere model and the kinetic Enskog-Landau equation for the charged hard sphere model in a 
pair collision approximation, respectively |4]. In the latter case, when <l>'(|r,j|) is the Coulomb potential 
of interaction, taking into account the features tq ^ -i-O, Aro ^ -i-O makes it possible to separate a colli- 
sion integral of the revised Enskog theory and a Landau-like collision integral in the limits t ^ -0 and 
T -oo, respectively. In the case of potential (8), in the region of tq +Q, Aro +G, o + Aro where 
the main contribution to a dynamics is defined by pair collisions of hard spheres, the memory function 
<Pnn(k;p,p'',t,t') can be calculated by expanding it over the density (a pair collision approximation), 
which was scrupulously done in papers by Mazenko ilit - fivll . 

Then, the kinetic equation Q) can be represented in the following form: 

|:/it(p;?) + -^/fc(p;0 = -^"/o(p)[c2(fc)-g2(fT)c2°(fc)] [ Ap'f'^{p';t) 
at m " m J " 

- j dp'tp^lit p, p')fi{p'; t) + innhih p)hfHt) 

t t 

-Jdp'l dt'e"'*-''W„[hp,p';t,t')fi{p';t')- j dfe'^'-''^ipnh(k;p;t,t']hfHt']. (9) 
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Here, 

><\foip'*)fk(P't)-foip')fiip*;t) + e'^'^-^''Mp'*)f'^{p'*;^ (10) 

is the Enskog-Boltzmann collision integral, where C2 (k) is the low-density limit of the direct correlation 
function and gzicr) is the pair distribution function. The step function 0_(x) is unity for x < and van- 
ishes otherwise. dD.a is the differential solid angle, 5 is unity vector. The precoUision and postcoUision 
momenta of the colliding hard spheres are denoted as (p, p') and {p* ,p'*), respectively. (p^n„{k; p, p'; t, t') 
is the part of the transport kernel related to the long-range interaction potential <l>'(|r;j|). Notably, the 
presented equation contains the Enskog-Boltzmann collision integral describing short-time dynamics of 
the hard sphere model. The collective effects related to the long-range interactions between particles are 
described by the functions iO„;,(fc; p), (p'„„(fc; p, p'; t, t'), (pnh(k) t, t') and by the equation for h^^^(t). Since 
the collective modes for the Enskog-Boltzmann model are well studied the investigation of time cor- 
relation functions and collective modes for the system of particles interacting through the potential (|8) 
turns out to be of great interest. In the case of the hard spheres system, the set of kinetic equations (2), l[9) 
reduces to the Enskog-Boltzmann kinetic equation Ulil . 

^/fc(p;f) + — /ifc(p;f) = -— «/o(p)[c2(fc)-g2(CT)c2(fc)] [ dp7fc(p';^) 

ot m m J 

-ng2{a)a^j dQ^ j dp' [d-\p-p']] 

X \Mp'*)hip; t) - h[p')fi[p*; t) + e^^-^'^/o(p'*)/j(p'*; f) - ^^-^^ MP) flip'; t)] . (ii) 

Projecting the Enskog-Boltzmann equation fiTt onto the first moments of the nonequilibrium one-particle 
distribution function a spectrum of collective excitations for the hard sphere model was obtained in 
[18]. Herewith, it is important to note that for the kinetic Enskog-Boltzmann equation we can consider two 
t5rpical limits: ka « 1 and fca » 1. In the hydrodynamic limit (ka « 1) the spectrum includes: heat mode 
zn{k) - -Dj^k^, where Dye is the thermal diffusivity coefficient in the Enskog transport theory |19]; two 
sound modes with eigenvalues given by z+ (fc) = +icfc - Fe fc^, where Fe is the sound damping coefficient 
and c is the sound velocity in the Enskog theory; two shear modes with eigenvalues given by Zy^ [k] - 
ik) = Zv(fc) = -v^k^, ve is the kinematic viscosity in the Enskog dense gas theory. In the limit ka » 1 
the Enskog-Boltzmann collision integral (Tot is transformed |11] into the Lorentz-Boltzmann collision 
integral which has only one eigenfunction ^i(p) - 1. Consequently, we obtain the diffusion mode only 
with the eigenvalue zj^lk) - -D^k^, where De is the self-diffusion coefficient as given by the Enskog 
dense gas theory. 

Let us now project the system of equations (2), (9) onto the first moments of the nonequilibrium 
one-particle distribution function. Thereafter, we perform simple transformations consisting in the tran- 
sition from the set of equations l|4) for averages ?z^(z), Jj(z),/j|'"(z),/i™'(-z) to the equations of general- 
ized hydrodynamics for averages bj,(z] - [«^(z),/j(z), h^iz) - h^'^(z] + h^^^{z)]. This permits to correctly 
define (see below) the generahzed viscosity coefficient via the transport kernel (13) and the heat con- 
ductivity coefficient via the transport kernel U.hfi(k,z). The averages ^j(z) satisfy the set of equations 
zbj_{z) - l.G{k; z)bi{z) - -{b^Xt - 0))'. In the limit ka » 1, the latter reduces to a single equation of 
diffusion for nj(z) in which the transport kernel "LgUc; z) - i\(p]:i^ (k)\^ i) corresponds to the Lorentz- 
Boltzmann collision integral flOt . In the opposite case, when ka « 1, the matrix t.Q(k;z) is defined as 
follows: Icihz) - liiihz) - iCluik) - nH(fc;z), 



tH{k;z): 



f in„j 



in;,^-n,,^ -{W3\cp^%m)-n 



hh 



(12) 

Ik.z) 
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Here we use the notations 

i:jj{k,z] = n,-,(fc,z)-2f^'(fc,z)[2^^"''^"(fc,z)]"'{in^;>(fc) + n^;f''"^(fc,z) (is) 
n,,H{k, z) = ^."•'^"(fc, z) + n^,"''"'(fc, z) + n;;^f "(fc, z) + uf/^Hk, z), (i4) 

where Z^^'^"(fc,z) = z - <^'^"(fc,z), X'^:''"^^,^) = n^^^'^-^fc-z) + n'^f "(fc,z) [l^;^'^"(fc,z)p 

xn^^'^"'(fc,z), 2'^"j(A:,z) = iO'^"j(A:)+n)^"^''^"(fc,z) [z^^''^"(fc,z)] i0^j>(A:). We can separate real and imag- 
inary parts in memory functions dlS) and lO as follows: Z ,• ,• (fc, z) = Z'. . (A:, &») + iZ'| . (fc, w) and Uhh {k, z) - 
U'^^ikyW) + in"^^(A;,fc)). Herewith, the contributions from the hard sphere dynamics with typical spatial- 
temporal scale To +0, Aro +0 are separated in the transport kernel (pnn(k; p, p'; t, t') only in the first 
term in the right-hand side of elements and hence in (13) . After these transformations we can ob- 
tain a spectrum of collective excitations in the hydrodynamic limit fccr « 1: heat mode ZH(fc) - -Djk^, 
where Dj is the thermal diffusivity coefficient for the system with the potential of interaction (8). It 
has the following structure: Dj - _Dte + Dj, Dj is determined through the corresponding elements Q 

of matrix of transport kernels (6). Dj - , Cp is a heat capacity at constant pressure, A' is the heat 

conductivity coefficient in the hydrodynamic limit: A' - limfc^o.w— o-^'C^.'^)- A'(fc,(y) is the generalized 
heat conductivity coefficient defined via elements of the matrix (6): A'(fc,w) - -^^-^n'^^Cfc.w), where 

cv{k) is the generalized heat capacity at constant volume dependent on the wave vector k; two sound 
modes z+(fc) = +ick- Tk^, where T is the sound damping and c - i-^p^siS)) sound velocity in 

the system with the potential of interaction (8), S(0) - S(k- 0), S(fc) is a static structure factor of the 
system with potential (8). r = j[Cplcv - 1)-Dt + jV^, where cy - cy^k = 0), 77^ - [^rj-^ + ri^) Imn is the 
longitudinal viscosity defined via the bulk viscosity rf' - rf^ + r]^ and the shear viscosity t?"*" = ?7j!r + t]^ 
coefficients. 77^ is the shear viscosity in Enskog theory, and r]^ is calculated in the hydrodynamic limit 
T]^ - limj;„o,ai— o'?['"(fc.'w)- '7['"(fc>^t') is the generalized shear viscosity coefficient defined via elements 
of the matrix (6) ?7T^(fc,w) = ^^I.".j-{k,ci)). Zf.(A:,z) is the transverse component of the generalized 

transport kernel l.jj{k,z), where the wave vector k is directed along the OZ axis. The longitudinal vis- 
cosity coefficient t]^ is calculated in the hydrodynamic limit r]^ - limj;^o,(u— o'7i'(fc.'^)> where r]^{k,a)) is 
the generalized longitudinal viscosity coefficient defined via longitudinal components of the generalized 
transport kernel Yjj{k,z): T]^[k,w) - ^^I.'j^(k,cD]; two shear modes with the eigenvalues given by 

Zvik) - -vk^. V = ve + vi is the kinematic viscosity v - rj^lnm for the system with the potential of inter- 
action (8). Here, v; is a contribution determined by the corresponding elements of the matrix of trans- 
port kernels l[6). In the limit fccr » 1, we obtain a diffusion mode, with the eigenvalue z^ik) - -D^k^, 
which is the same as in the Enskog theory. 

As we can see from the above expressions, presence of the long-range part in the potential of in- 
teraction entails a renormalization of all the damping coefficients in the collective modes spectrum. In 
particular, contributions related to long-range potential appear in heat and sound modes as well as in 
shear modes. Nevertheless, diffusion mode remains unchanged. 



4. Conclusions 

In this brief report within the framework of consistent description of kinetic and hydrodynamic pro- 
cesses we considered a set of kinetic equations for the potential of interaction of the system presented by 
the sum of hard spheres potential <l'''*(|f,y|) and a certain smooth one <l''(l?;jl)- In this case, we separated 
the Enskog-Boltzmann collision integral describing a collision dynamics at short distances from the colli- 
sion integral of the kinetic equation for the nonequilibrium distribution function. Applying the procedure 
of projecting onto the moments of the nonequUibrium distribution function to the equations 0, l[9) we 
obtain a set of equations for hydrodynamic variables. Based on this set of equations a spectrum of collec- 
tive excitations was obtained in the Umits ka « 1 and fccr » 1. We showed that, besides the contribution 
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from the hard spheres potential, all hydrodynamic modes contain contributions from the long-range part 
of potential. These contributions make the damping coefficients closer to the ones known from the hy- 
drodynamic theory. Here, we formally presented the contribution from the long-ranged part of potential, 
since the latter, for example the Coulomb one, will contribute into the transport kernels l[3). Moreover, we 
can separate the linearized Landau-like collision integral describing pair collisions in (pnnik; p,p'; t, t'), 
whUe (phnik; p; t, t'), (pnhik; p; t, t'), (phhik; p; t, t') take into account collective Coulombic interactions. 
Evidently, calculation of the elements l[7) of matrix tl{k; z) wiU depend on the model of time dependence 
(exponential, Gaussian etc.) for transport kernels (3). When a spectrum of collective excitations is known, 
a whole set of time correlation functions can be investigated. In particular, it makes possible to investigate 
the behaviour of the dynamic structure factor and, in the case of potential (8), to separate a contributions 
from the hard spheres potential and the long-range part of potential in it. 
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f\o npo6yieMi/i y3roA>KeHoro oni/icy KiHeri/mHi/ix ra 
riApoAi/iHaMiMHi/ix npoi^ecie y rycTi/ix ra3ax Ta piAi/iHax: 

CneKTp K07ieKTI/IBHI/IX 36yA>KeHb 

B.B. MapKiB, I.n. Ouen^H, M.B. ToKapnyK 

iHCTwryr (|)i3HKW KOHfleHCOBaHwx cucreM HAH YKpaiHH, sy/i. L CBeHi^ii^bKoro, 1, 7901 1 /IbBis, YKpaiHa 

Ha OCHOBi OTpWMaHWX paHilUe ySara/lbHeHHX HeMapKiBCbKWX piBHSHb fl/ia HepiBHOBaXHOI OflHO^aCTUHKO- 

BOT (|)yHKL^iT po3nofli;iy ra cepeflHboro SHaneHHA rycTWHW noTeHi^ia/ibHoT 4acTi/iHH eHra/ibnii [Markiv B.B., 
Omelyan LP, Tokarchuk M.V., Condens. Matter Phys., 2010, 13, 23005] floc/iiflxyerbcfl cneKip KO/ieKTHBHUx 36y- 
flxeHb, KO/in noTeHL4ia;i BsaeMOfliT MixnacTUHKaMi/i npeflcraB/ieHO cyMOio noTeHi4ia;iy TBepflux cefiep la fleaKoro 
fla/ieKOCflXHOro noreHL^ia/iy. 

K/iK>HOBi cjioBa: KiHertiKa, riflpoflUHaMiKa, KineTi/mHi pibhrhhr, <pyHKu,ii naM'RTi, KoneKTi/iBHi MOfli/i 
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